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1 Introduction 

Recently several factors converged to ignite a renaissance in charm physics. These 
, factors include: 

O \ • the need for precision measurements from the charm sector to interpret quantita- 

^ ■ tively CP violation results from the beauty sector, 

' • technical and computational developments in Lattice QCD that led to precise 

o ■ 

■ • discovery of new - generally unanticipated - charm meson states {e.g., D*q(2317) 

Oh: 



calculations that can face precision experimental challenges in the charm sector, 
discovery of new - ge 
and D,i(2460)), and 



^ ■ • substantial increases in the precision and reach of charm decay measurements due 

to production of much larger charm data sets (particularly at CESR, KEKB, and 
PEP II) and excellent detectors in the experiments at these facilities. 
^ \ In this report I describe recent results in hadronic decays of the stable D mesons 

' and -Ds) and Dalitz analyses of resonance structure in their hadronic decays. 

These results come from experiments at electron-positron colliders (BaBar, Belle, 
and CLEO-c) and fixed target experiments (E791 and FOCUS). 



2 Absolute D and Branching Fractions 

CLEO is providing new precise measurements of absolute and branching 
fractions using and decays from e"''e~ — ?/'(3770) D^D^ or D^D . The 
mass of the '0(3770) is below the threshold for DDtt decays, so no additional pions are 
produced. Furthermore, the multiplicities of and decays are low, so events are 
extremely clean. CLEO measures leptonic, semileptonic, and key hadronic branching 
fractions using a double tagging technique pioneered by MARK III [1, 2]. Most 
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other D branching fractions [3] are measured relative to a reference mode, usually 
K~7T^ or K^n^TT^. CLEO has published absolute branching fractions 

for key Cabibbo-Favored Decays (CFD) to hadrons obtained from 56 pb"^of data 
[4]. I report here a preliminary update from 281 pb~^ of data; this is the first public 
presentation of these results. Some other branching ratios utilizing the 281 pb~^ data 
sample have been published or submitted for publication. 

The MARK III double tag technique derives absolute branching fractions from 
measurements of Single Tag (ST) and Double Tag (DT) yields. In ST events, the 
D{D) is reconstructed in a specific final state, while the decay products of the D{D) 
are not observed, e.g., D ^ i and D —>■ X. In DT events, both the D and the D 
are reconstructed in specific final states, e.g., D ^ i and D ^ j. ST and DT yields 
are given by A^'j = A''^^? = ^dd ^1^1^ ^ij ~ ^dd respectively, 

where: i and j are the decay modes of the D and D, respectively; Bi, and are the 
corresponding branching fractions; e^, ej, and e^j are the ST and DT efficiencies for 
these modes; Ni, Aj, and Ajj are the corresponding ST and DT yields; and A"^-^ is the 
number of D^D or D^D~ events produced in the experiment. Branching fractions 
and Njj-^ can then be obtained from 

where ej is the efficiency for detecting D ^ j ST events, and CP symmetry has been 
assumed so = Bi. Obviously the luminosity is not required to determine Bi or 
Njfij, but the cross section for e~^e~ — > '?/'(3770) — DD can be determined from A^^-^ 
and the luminosity. Furthermore, ~ e^ej so the branching fractions and number 
of events obtained by this method are relatively insensitive to uncertainties in the 
efficiencies. 

CLEO utilizes a fit [5] that obtains all D^ and D^ branching fractions, as well 
as the numbers of D^D^ and D^D~ events from a simultaneous fit to all ST and DT 
D^ and D^ yields. All statistical and systematic uncertainties and their correlations 
are properly taken into account in this fit. 

Candidate D or D mesons are selected using mode-dependent requirements on 
AE = E{D)—Eo, the difference between the measured energies E{D) of the candidate 
and the beam energy Eq. Then ST and DT yields are obtained from fits to beam 
constrained mass Mbc distributions, where the beam energy Eq is substituted for the 
measured energy E{D) oi a. D or D candidate. Substitution of Eq for E[D) improves 
substantially the mass resolution of the D candidates. The Mbc distributions for 
the DT events are illustrated in Figure 1. Clearly these Mbc distributions are very 
clean with little background. Summed over all modes, the yields obtained from the 
CLEO-c 281 pb-^ sample are 230,225 ST and 13, 575 ± 120 DT D^'^ events, and 
167,086 ST and 8, 867 ± 97 DT D^D~ events. With these yields, systematic errors 
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dominate; to be conservative for these preliminary results, most systematic errors are 
essentially the same as those determined for the 56 pb~^ results. Intensive study of the 
systematic uncertainties are underway and CLEO expects to be able to substantially 
improve most of them. 
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Figure 1: Histograms of Mbc distributions of D and D candidates in double-tag 
events from 281 pb'^ of CLEO-c ?^(3770) data; (left) D+D~ DT candidates and 
(right) D^D candidates. 



Mode 


B (%) 


D^^ 


K-71+ 




3.87 ±0.04 ±0.08 


D^^ 


K-TT + 




14.6 ±0.1 ±0.4 


D^^ 






8.3 ±0.1 ±0.3 


D+ ^ 




-71 + 


9.2 ±0.1 ±0.2 


D+ ^ 






6.0 ±0.1 ±0.2 


D+ ^ 






1.55 ±0.02 ±0.05 


D+ ^ 






7.2 ±0.1 ±0.3 


D+ 






3.13 ±0.05 ±0.14 


D+ ^ 


K+K 


-71 + 


0.93 ±0.02 ±0.03 



Table 1: Preliminary CLEO results for hadronic and D+ branching fractions, with 
their statistical and systematic errors. These results come from 281 pb~^ of CLEO-c 
■0(3770) data, with conservative systematic errors. 
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The preliminary CLEO results for three and six D"*" hadronic decay modes are 
given in Table 1 and compared to PDG04 [6] averages in Figure 2. It is clear that the 
preliminary CLEO results are substantial improvements over previous measurements. 
PDG04 averages were used for the comparison because CLEO-c 56 pb~^ results were 
included in the PDG06 averages. Final State Radiation (FSR) is included in the 
CLEO-c Monte Carlo (MC) simulations for determining efficiencies. If FSR had not 
been included in the simulations the branching fractions would decrease by ^ 2%. 
Most other measurements of hadronic D branching fractions did not take FSR into 
account. 
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Figure 2: Ratio i3(CLEO)/i3(PDG04) of preliminary CLEO-c hadronic D° and D+ 
branching fractions to the 2004 PDG averages. The widths of the PDG bars cor- 
respond to the errors in those averages. The CLEO-c points have statistical and 
systematic error bars. 



3 Cabibbo Suppressed and Decays 

3.1 Singly-Cabibbo-Suppressed and D'^ Decays to Pions 

Using the full 281 pb~^ 11^(3770) data sample, CLEO measured Singly-Cabibbo- 
Suppressed Decays (SCSD) of and to multipion final states [7]. The Mbc 
distributions are illustrated in Figure 3. The branching fractions in Table 2 were 
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Figure 3: Histograms of Mbc for (plots on the left) of {n)7i events and (plots on 

the right) {n)^ from the 281 pb~^ '?/'(3770) CLEO-c data sample. Histograms 

of Mbc for the reference modes K~'k^ and K^-k^-k^ are included with 

the D n{7i) histograms on the left and right, respectively. 



Mode 


CLEO-C B (10-3) 


PDG04 B (10-3) 


TT + TT 




1.39 ±0.04 ±0.03 


1.38 ±0.05 






0.79 ±0.05 ±0.04 


0.84 ±0.22 






13.2 ±0.2 ±0.5 


11 ±4 




TT 


7.3 ±0.1 ±0.3 


7.3 ±0.5 






9.9 ±0.6 ±0.7 








4.1 ±0.5 ±0.2 








1.25 ±0.06 ±0.08 


1.33 ±0.22 






3.35 ±0.10 ±0.20 


3.1 ±0.4 






4.8 ±0.3 ±0.4 








11.6 ±0.4 ±0.7 








1.60 ±0.18 ±0.17 


1.82 ±0.25 



Table 2: CLEO branching fractions for and decays to multiple pions. 



obtained from measurements of i3(i5° — n{7i))/B{D'^ — > K~7r~^) and B{D^ — >■ 
n{7r))/B{D~^ K~7T~^7r~^). The reference branching fractions used to determine 
the reported multipion branching fractions were B{D'^ K~n~^) = (3.84 ± 0.07)% 
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and B{D^ K~tc'^tc^) = (9.4 ± 0.3)%. These were obtained by averaging the 
56 pb~^ CLEO-c results and the PDG04 averages. These new CLEO measurements 
are compared to the PDG04 [6] averages in Figure 4. 

■ PDG04 ■ CLEO-c 281 pb"' 0140107-003 

Jl+71" 



7l*7t*7t 71 K* 

0.5 1 1.5 2 

'B(CLEO)/'B(PDG04) 

Figure 4: Ratio i3(CLE0) /i3(PDG06) of CLEO-c measurements of D ^ 71(11) branch- 
ing fractions to the 2006 PDG averages. The widths of the PDG bars correspond to 
the errors in those averages. 

As described in Section 3.4, BaBar has also reported a measurement of B{D^ 
TT+TT^) along with the first observation of the Doubly-Cabibbo-Suppressed-Decay 
(DCSD) branching fraction B{D^ K'^-k'^). The BaBar result for the ratio of 
B{D^ — >■ tt+tt") to the reference branching fraction B{D^ — s> K~tt^ti^) is in excellent 
agreement with the CLEO measurement. 

The ratio of the isospin amplitudes Aq and A2 for D ^ nir decay and their phase 
difference can be determined from the values of Bln^ir"), BItt^tt^), and B{7r~^'K^) and 
the and decay widths [8]. Only 7 = and 2 states are allowed for a two pion 
system, since an isovector state is forbidden. Using D lifetimes from PDG04 [6] and 
the branching fractions in Table 2, CLEO finds A2/A0 = 0.420 ± 0.014 ± 0.001 and 
the relative phase S = (86.4 ± 2.8 ± 3.3)°. In contrast to the isospin amplitudes in 
K ^ TTir decays, the two isospin amplitudes in D ^ tttt decay are comparable [8]. 
Furthermore, the large phase difference between the two amplitudes indicates that 
final state interactions are important. 
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3.1.1 Searches for 77 and uo in multipion and decays 



Untangling the resonant substructure of these decays with more than two pions would 
require systematic Dalitz analyses. However, searches for 77s and ujs among the decay 
pions is a first step in this direction. The technique used is illustrated in Figure 5, 
where three-pion invariant masses, M(7r+7r~7r°) are plotted for D 01 D candidates 
within the original AE requirement and in sidebands. Although combinatorial back- 
grounds are large, there are significant signals above background peaks in: 777?° in 
7r^7r~7r°7r°, r/vr"'' in tx^ ii^ tt^ , and uj'k^h~ in 7r+7r~7r'^7r~7r"^. 

The branching fractions for these three channels and 90% UL confidence intervals for 
the other three modes are given in Table 3. The branching fractions for ?77r"'" and 
uoTi'^Ti" are significant fractions of the total branching fractions for the parent decays. 
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Figure 5: M(7r+7r ir^) distributions for multipion and decays from CLEO-c 
data. 



3.2 Singly-Cabibbo-Suppressed Decays from BaBar and Belle 

BaBar [11] and Belle [12] have used their enormous data samples to measure pre- 
cisely the ratio B{D'^ 7i^7i~7i^)/B{D^ — >• K'lr^n^) of a singly-Cabibbo-suppressed 
hadronic decay to a similar Cabibbo-favored hadronic decay. The BaBar and Belle 
data are illustrated in Figures 6 and 7, respectively. The invariant masses of D can- 
didates are used in these analyses, since these data were taken at the T{4S) where 
many DD channels with multipions are open (rather than only the D^D^ and D^D~ 
channels at the ip{3770)) so beam-constrained masses cannot be used. Therefore the 
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Mode 


B (10-3) 


PDG04 B (10-3) 


?77r° 


0.62 ±0.14 ±0.05 


— 


rjTi^ 


3.61 ±0.25 ±0.26 


3.0 ±0.6 




< 1.9 (90% CL) 






< 0.26 (90% CL) 






< 0.34 (90% CL) 






1.7±0.5±0.2 





Table 3: Branching fractions and upper limits for rj and uo production in multipion 
and decays from CLEO-c data. 




Figure 6: BaBar data for (left) — > K tt+tt^ decays and (right) ir^ir n' 

decays. [9, 10] 



mass resolutions are substantially worse in the BaBar and Belle data, than in the cor- 
responding CLEO-c data (compare Figures 6 and 7 with Figure 5. However, signal 
to background ratios in the BaBar and Belle data are substantially better than the 
corresponding ratio in the CLEO-c data. The results of the three measurements are 
given in Table 4. Some of the strengths and weaknesses of these measurements are 
apparent from the yields, errors, and luminosities in the table. Note that the Cabibbo 
suppression of BlD'^ — > tt+tt-tt'') is about a factor of 10 relative to B{D^ K~7t'^-k^), 
which is the hadronic decay mode with the largest branching fraction. 

In the same analysis BaBar measured B{D^ K~ K^tx^)/B{D^ K^tt^ti^). 
The data are illustrated in Figure 8 and the result is included in Table 4. The BaBar 
result is significantly larger and much more precise than a previous CLEO result from 
a much smaller data sample [13]. 
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Figure 7: Belle data for (left) Ji^Tr+Tr^ decays and (right) — > 7r'^7r~7r' 

decays. 





Yield 


B{D'> 7r+7r-7rO)/i3(DO ^ J^-vt+ttO) 


Luminosity 


BaBar 


60, 426 ± 343 




^10.59 ±0.06 ±0.13) X 10^2 


232 fb-i 


Belle 


22, 803 ±203 




( 9.71 ±0.09 ±0.30) X 10-2 


357 fb-i 


CLEO-c 


10, 834 ± 164 




( 9.01 ±0.18 ±0.39) X 10-2 


281 pb-i 








^ K-K+n^)/B{D^ A'-tt+ttO) 




BaBar 


10, 773 ± 122 




(2.37 ±0.03 ±0.04) x lO'^ 


232 fb-i 


CLEO II 


151 ± 42 




(0.95 ±0.26) X 10-2 


2.7 fb-i 



Table 4: Measurements of SCSDs from BaBar, Belle, and CLEO. The Belle result is 
preliminary, while the BaBar and CLEO results have been published. The CLEO-c 
branching ratio is calculated from measured ratios of both B{D^ -k^-k^-k^) and 
B{D^ — > K'n'^ii^) to B{D^ — > K~7r~^), without taking correlations into account. 
Proper use of correlations would reduce the systematic error. 



3.3 Multi-Kaon Modes from FOCUS 

FOCUS studied decays to final states with two or three charged or neutral 
kaons [15, 16]. The data for the reference mode — >■ A'^Tr+Tr" and two of the 
neutral kaon signal modes are illustrated in Figure 9 and the ratios of signal branch- 
ing fractions to the reference branching fraction are given in Table 5. This is the first 
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m (K"K+JiO) [GeV/c^] 



Figure 8: BaBar data for D° R-R+ti^ decays. [9, 10] 




Figure 9: FOCUS data for (left) ^ Ji^vr+Tr- decays, (middle) ^ R^R^R^n^ 
decays, and (right) R^^R^gR^g decays [14]. 



observation of the modes R^R^R'^tt^ , which are combined to one branch- 

ing fraction. The table also includes a measurement of the branching ratio B{D^ — *• 
R~ R~^7T~^7r~) / B{D^ —>■ R~7T^7r^7T~). This branching ratio was measured in a sep- 
arate FOCUS analysis of resonant substructure in the the four-body decay D'^ 
R~ R'^TT^TT" described in Section 6.3. 

From the table it is evident that Cabibbo-suppression, relative to the Cabibbo- 
favored reference mode, is on the order of a factor of 100 for the two-kaon decay modes. 
However, the Cabibbo-favored decays to three kaons are suppressed by comparable 
factors. Presumably this is due to the necessity of popping an ss quark-antiquark 
pair to produce the additional two kaons. 
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Mode 


FOCUS i3^,ode/^ref (%) 


D^^ 




1.44 ±0.32 ±0.16 


SCSD 


D^^ 


K'^KItt+tt- 


2.08 ±0.35 ±0.21 


SCSD 


D^^ 




1.06 ±0.19 ±0.10 


CF 




KIKIKI 


1.79 ±0.27 ±0.26 


CF 




K-K+TT+TT- 


2.95 ±0.11 ±0.08 


SCSD 



Table 5: FOCUS branching fractions for Cabibbo favored (CF) decays and singly- 
Cabibbo suppressed D'^ decays to multiple kaons. For the neutral kaon modes, the 
reference branching fraction is B{D^ —>■ K^ti^ti'), while the reference branching 
fraction for the charged kaon mode is B{D^ K^ir^Ti^n'). 

3.4 The Doubly-Cabibbo-Suppressed Decay — > K'^7r~ 




BaBar has reported the first observation of the DCSD K^tt^ [18] and the 

precise measurement of a branching ratio for the SCSD decay tt^tt^ mentioned 

in Section 3.1. The BaBar data are illustrated in Figure 10, and the CLEO data [19] 
for the same modes are illustrated in Figure 11. Note that the mass range of the 
CLEO data is much narrower than the mass range of the BaBar data, illustrating the 
improvements in mass resolution gained by being able to use beam-constrained masses 
instead of invariant mass. The ratios of the signal branching fractions B{D^ tc^tt^) 
and B{D^ K'^ii^) of the signal modes to reference branching fraction B{D^ 
K'TT^TT^) are given in Table 6. BaBar's ability to use tighter cuts on their enormous 
data sample to compensate for the better mass resolution of the CLEO-c data is 
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Figure 11: CLEO-c data for (left) D+ K-TT+n+ decays, (middle) D+ vr+vr' 
decays, and (right) Z)+ K~^7r~ decays. 



evident from the figures and results. The Cabibbo suppression of the K~^it' 
mode is close to a factor of 5. 





D+ tt+ttO (SCSD) 


D+ ^ K+n^ (DCSD) 




^mode/^ref 


^mode/^ref 


BaBar 


(1.33 ±0.11 ±0.09) X 10-2 


(2.68 ±0.50 ±0.26) x lO^^ 


CLEO-c 


(1.33 ±0.07 ±0.06) X 10-2 


(2.40 ±0.38 ±0.16) X 10-^ 



Table 6: The ratios of the branching factions for the SCSD decay tt^tt^ and 

the DCSD decay K^ir^ from BaBar and CLEO-c, to the branching fraction 

for the reference mode D'^ K^tt^tt^. The CLEO-c result for tt+vt'^ is from 

the SCSD analysis in Section 3.1. 



3.5 Comparison of D — > ^s'^ ^ — ^ ^l'^ Decay Rates 

Cabibbo-Favored and Doubly-Cabibbo-Suppressed amplitudes contribute to the de- 
cay D — > K^TT. The observed final states are D KgTi and D Kj^tt, and Bigi and 
Yamamoto pointed out that interference between CF and DCS amplitudes can lead 
to different rates for D Kgir and D K^tt [20]. CLEO has measured prelimi- 
nary branching fractions for these decays by fully reconstructing D — s> Kgir decays 
and reconstructing D K^ir decays using missing masses [21]. The CLEO data 
for D — > Kill decays are illustrated in Figure 12. The preliminary results given in 
Table 7 are reported in terms of the ratio R{D) defined by: 

^ B{D ^ KgTT) - B{D ^ KItx) 
^ '~ B{D^ K^Tr) ± B{D ^ Kin) 
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The measured value of R{D^) is consistent with zero, while R{D^) is significantly 
larger than zero. U-spin and SU(3) predict [22] R{D°) = 2tan2(6'J giving R{D°) = 
0.109 ±0.001, which is in good agreement with the experimental result. On the other 
hand, R{D~^) is not so simple because internal spectator diagrams contribute to both 
— >■ K^7r~^ and K^tt^ , but external diagrams contribute to the former and 

annihilation diagrams contribute to the latter. 




Figure 12: CLEO data for (left) D'^K^ti^ decays and (right) D+ KItx+ decays. 



R{D^) 
RiD"^) 



0.030 ±0.023 ±0.025 
0.122 ±0.024 ±0.030 



Table 7: CLEO-c Preliminary measurements of the normalized branching fraction 
differences between D Kl'^ decays and D — ^ Ks"^ decays. 



4 Absolute Dg Branching Fractions 

Previously, absolute Dg branching fractions have not been well known, primarily due 
to difficulties of determining the total number of Dg mesons produced in the experi- 
ment [23, 24]. CLEO has now addressed this problem with a double tag measurement 
of absolute Ds branching fractions in e~^e~ collisions just above the DfD*"^ thresh- 
old [25]. 

4.1 The Ds Production Cross Section 

Although measurements of the e+e~ annihilation cross section, a{e^e~), at energies 
above Ecm = 3.8 GeV have existed for some time [3], little was known about the 
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composition of the final states. In order to find a favorable point for studying DT 
DfD^ events, CLEO scanned that region with integrated luminosities of ~ 5 pb"^ per 
point and fast turnaround of results for feedback. After the scan, more luminosity 
was accumulated at or near Ecm = 4.17 GeV where the cross section for DfD*^ 
production peaks with a{e'^e~ DfDl"^) 0.9 nb. Preliminary measurements [25] 
of cross sections for producing various DD, DD , and D*D pairs are illustrated in 
Figure 13. CLEO used a total of 195 pb~^ of data at Ecm = 4.17 GeV in the analysis 
of Ds decays described below. 




4.1 4.15 

(GeV) 



4.1 4.15 

(GeV) 



Figure 13: Preliminary CLEO-c cross sections for production of (left) D^*)!)''*'' and 
(right) Df^*^D~^*^ above the threshold for DD production in e~^e~ annihilation. In 
each graph thresholds are indicated by the vertical dashed lines. 



4.2 Analysis of DfDf^ Events 

CLEO chose to ignore the 7 or vr'' from D* decay to avoid the low efficiency for 
detecting the low energy 7 or 7r° and the uncertainties in the efficiency for detecting 
either one of them. Instead, CLEO selects DfD*"^ events using the invariant mass 
M[Ds) of the Dg candidates and their beam-constrained mass Mbci which is actually 
a proxy for momentum of the candidate. The Mbc distribution for the Dg candidates 
that were produced directly in the annihilation is relatively narrow, while the M^c 
distribution for D^ candidates resulting from decay is much broader. This is 
illustrated in Figure 14. Neither distribution is centered exactly on the Ds mass, 
because the beam energy is not exactly the energy of either the direct Dg or the 
average energy of the Dg from Z)* decay. 
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Figure 14: Histogram of Mbc for Df — >■ K~K^7r^ events from CLEO-c data. The 
narrow peak at the Mbc = 2.04 GeV is from Df produced directly, while the broad 
peak from Mbc ~ 2.01 GeV to Mbc ~ 2.07 GeV is from Df produced in D*^ decay. 



CLEO determines single-tag yields by fitting M{Ds) distributions for candidates 
that pass a very loose cut on Mbc- Double-tag yields are determined by counting the 
numbers of events in signal regions in the M{D~) vs. M{D^) plane, and subtracting 
backgrounds determined from the numbers of events in sideband regions. These 
procedures are illustrated for Df K~K^ti^ candidates in Figure 15. 

The fit [5] procedure used in the measurements of absolute and branch- 
ing fractions is applied to the ST and DT yields for six Dg decay modes. Preliminary 
absolute branching fractions obtained from this analysis are given in Table 8, and 
these branching fractions are compared to the PDG06 averages in Figure 16. The 
preliminary CLEO-c results are clearly significantly more precise than the PDG06 
averages. CLEO has an additional 130 pb~^ of data at this energy to be analyzed. 
These data will be included in the publication of the 195 pb~^ data sample. 

Belle measures B{D^ K~ K^tt^) by applying a clever partial reconstruction 
technique to e+e~ — > 0*^^ D^i{2b?>Q) events [26] that produces double tag data sam- 
ples. A total of 552.3 fb~^ of data taken at the T(4S') were used in this analysis. 
In the first tag, the D^i is fully reconstructed in the Dsi — > D*K mode, while only 
the 7 from D* — ^ Dg'-y is observed. In the second tag, the is fully reconstructed 
in the decay chain D*^ Df'-j followed by Df — ^ K~K^7r^, while only the kaon 
from Dsi D*K is observed. Signal yields are determined by fitting missing mass 
distributions for the particles that are not reconstructed. The yields from the first tag 
are proportional to the branching fractions for the D* decay chain, while the yields 
from the second tag are proportional to the branching fraction for Df — >■ K~ K^tt^. 
Hence, B{Df — > K~K~^7r~^) can be obtained from the ratio of the second tag yield 
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M(Df) (GeV/c^) M(D:) (GeV/c') M(Dp - M(D;) (GeV/c') 

Figure 15: Plots illustrating the CLEO-c Df decay data: (left) the invariant mass 
distribution for single-tag Df — >■ K~K~^7r'^ candidates, (middle) two-dimensional his- 
togram of M(D7) vs. M(D+) for double-tag Df R-R+ti^ candidates, (right) the 
mass difference M[D^) — M{D~) for double-tag Df K~ K^-k"^ candidates. In the 
two-dimensional histogram, the narrow (blue) box around the mass peak illustrates 
the signal region, and the two wider (red) boxes further away from the peak illustrate 
the sideband region. 



Dj Mode 


B (%) 


KsK+ 


1.50 ±0.09 ±0.05 


K-K+Tr+ 


5.57 ±0.30 ±0.19 


K-K+Tx+n^ 


5.62 ±0.33 ±0.51 


7r+7r+7r" 


1.12 ±0.08 ±0.05 




1.47 ±0.12 ±0.14 




4.02 ±0.27 ±0.30 



Table 8: Preliminary absolute branching fractions for six hadronic Df decay modes 
measured by CLEO-c. 



to the first tag yield, using known branching fractions for the D* decay chain, and 
efficiencies. 

The preliminary Belle and CLEO-c results are compared in Table 9 to each other 
and to the PDG06 [3] average. With its large error, the PDG 06 value is consistent 
with the Belle and CLEO-c results. However, the more precise Belle and CLEO-c 
are not very consistent with each other. Hence, more work will be required to be 
confident that we know absolute Df branching fractions. 
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PDG 2006 fit BCLEO-c 195 pb ' Preliminary 




' ' ' I ' ' ' I ' ' ' ^ ' ' ' I ' ' ' 
0.4 0.6 0.8 1.0 1.2 1.4 

'B(CLEO)/'B(PDG06) 

Figure 16: The ratios i3(CLEO)/i3(PDG06) of preliminary CLEO-c hadronic 
branching fractions to the PDG06 averages. The widths of the PDG bars correspond 
to the errors in the averages. 





B{Dt ^ K-K+'K+) (%) 


Belle 


Preliminary 


4.0 ±0.4 ±0.4 


CLEO-c 


Preliminary 


5.57 ±0.30 ±0.19 


PDG06 




5.2 ±0.9 



Table 9: Comparison of preliminary Belle and CLEO-c measurements of the absolute 
branching fraction for K~K~^7r~^ decay with the PDG06 average. 



4.3 Partial Df K K^tt'^ Branching Fractions 

The branching fraction of the decay chain — » 07r+ — > K~K~^7r'^ is one of the 
largest branching fractions, and - indeed - it was the decay sequence that provided 
the first observation of the Df. A branching fraction called B{Df — > ^vr"*") has often 
been used as the reference branching fraction for Ds decays. As in the discovery, 
it is usually derived from a narrow mass cut around the (p peak in the M[K~^K~) 
distribution in K'K'^n^ events. 
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However, E687 [27] has published and FOCUS [28] has reported significant con- 
tributions from /o(980) (or ao(980)) in the (pn region of the K^K'^tt'^ Dahtz 
plot. Figure 17 illustrates the invariant K'^K~ mass distribution for D'^ K~K~^7i^ 
events. A noticeable contribution (presumably scalar) is visible under the (p peak in 
in this distribution. These contributions from other processes constitute approxi- 
mately 5% of the total yield in a reasonable mass range in M{K~ K^) around the 
(j) mass, M^. These contributions, which are not from D'^ (pTi^ decays, are com- 
parable to current CLEO-c errors for partial branching fractions obtained using a 
narrow mass cut around M^. Hence, it is now difficult to make a case for iden- 
tifying such a partial branching fraction with B{D'^ — > 07r+). In response to this 
dilemma, CLEO now reports partial branching fractions. Bam = B{D^ K~K'^tt'^) 
with \M{K~K~^) — M^\ < AM MeV/c^. Preliminary CLEO-c measurements with 
AM = 10 and 20 MeV/c^ are given in Table 10. I advocate using partial branch- 
ing fractions like these for the reference branching fractions for other decays. 
As D'^ branching fraction measurements become more precise, I hope that we can 
reach agreement on a reasonable choice for a Df reference branching fraction for 
most other decays that is not called B{D^ 07r+), even though it has non-negligible 
contributions from other resonances in the K~K~^7t~^ Dalitz plot. 




K+K- mass (GeV/c^) 



Figure 17: The invariant mass distribution M[K^K~ 
dates in CLEO-c data. 



for D: 



candi- 



5 Inclusive D , D , and Ds Decays to ss 

CLEO has measured branching fractions for inclusive -D", D^, and Dg decays to 
■qX, r]'X, and (j)X [29]. Due to the ss content of rj, t]', and </>, we expect that the 
inclusive branching fraction for Dg decays to these particles would be larger than the 
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Bam (%) 


Bio 


1.98 ±0.12 ±0.09 


B20 


2.25 ±0.13 ±0.12 


PDG06 


1.77 ±0.44 



Table 10: Preliminary partial branching fractions Bam = B{Df K K^ti^) with 
\M{K-K+) - M^\ < AM MeV/c^ measured by CLEO. The PDG06 resuh is deter- 
mined from B{(f) K~K^) and reported B{Df (p7i~^) values. 



branching fractions for corresponding and decays. CLEO fully reconstructs 
one D in a favorable mode used in the absolute measurements of hadronic branching 
fractions. Selection criteria for these D decay candidates include requirements on AE 
(see Section 2). Then CLEO searches for 77, 77' and (j) in the decay products from the 
other D. For this analysis, CLEO utilizes 281 pb"^ of ?^(3770) data for D° and D+ 
decays and 195 pb~^ of data at Ecm ~ 4.17 GeV for decays. 

0140107-004 




0.35 0.40 0.45 0.50 0.35 0.40 0.45 0.50 0.35 0.40 0.45 0.50 



M(Ti') - M{t|) (GeV) M(Ti') - M(ti) (GeV) M(ti') - M(ti) (GeV) 

Figure 18: Histograms of the mass differences between t]' and rj candidates with 
7]' — s> rjiT'^TT" with 7] 77: (a) candidates for — > r]'X decay, (b) candidates for 
JJ+ _^ j^'x decay, and (c) candidates for rj'X decay. 

Figure 18 illustrates the data used to determine the inclusive ri'X signal yields. 
For these inclusive modes, CLEO uses distributions in M{rj') — M(?7), the difference 
between the invariant masses of the rj' candidate and the rj candidate in its decay 
chain. Similar distributions for M{ri) and M{(f)) are used to determine the signal 
yields for r]X and (pX. Fitted yields from AE sidebands are then subtracted from 
the signal yields. The resulting inclusive branching fractions are given in Table 11. 
These results lead to several qualitative observations including: 
• r/'s and 0s are relatively rare in and decay. 
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• ?7S with their lower mass and larger light quark content are produced at substan- 
tially higher rates in and decays than 77' and 0, 

• the (f) rate in Dg decay is much higher than it is in and decay, and 

• the higher (p rates in Ds decay can be used to separate Dg from and at 
T(55') and hadron colliders. 

CLEO has already used the observation in the last bullet to study the Bg fraction in 
the final states in T{5S) decay [30]. 



Mode 


B{D') (%) 


B{D+) (%) 


B{Dt) (%) 


7]X 


9.5 ±0.4 ±0.8 


6.3 ±0.5 ±0.5 


23.5 ±3.1 ±2.0 


T]'X 


2.48 ±0.17 ±0.21 


1.04 ±0.16 ±0.09 


8.7 ±1.9 ±0.8 


(j)X 


1.05 ±0.08 ±0.07 


1.03 ±0.10 ±0.07 


16.1 ± 1.2 ± 1.1 



Table 11: Branching fractions for inclusive D ssX decays measured by CLEO. 



6 Dalitz Analyses of Hadronic D Decays 

In the last decade, large samples of hadronic D decays have been available to the 
E791 and FOCUS collaborations from fixed target experiments at Fermilab, and to 
the CLEO, BaBar, and Belle collaborations from e~^e~ collider experiments. Dalitz 
analyses of these data samples have enabled new insights into the resonant substruc- 
ture of multibody hadronic D decays. Historically, simply disentangling the resonant 
substructure has been the major emphasis of these studies. More recently Dalitz 
analyses of D decays have been stimulated by the possibility that the results can be 
used to determine the CKM matrix angles 7 or 03 . 

6.1 Dalitz Analyses of — ^ tt+tt+tt" Decays 

It has long been understood that several TT~^7r~ resonances are likely to contribute to 
the resonant substructure of D'^ tt+tt+tt" - decay. These include p°(770), /o(980), 
f2{g), and possible other higher resonances. On the other hand, there has been 
evidence (not always conclusive) that a lower mass S-wave 7r~^n~ state called the a'^ 
might exist and might contribute to — >■ 7r"'"7r"'"7r~ decay. Finding proof for or 
against the existence of this state and its contribution to — > vr+Tr+vr" decays has 
been a major concern of Dalitz analyses of this final state. 

The choice of the theoretical model used to fit the Dalitz plot has complicated 
interpretation of these analyses. The isobar model, in which the decay amplitude is 
written as a sum of Breit-Wigner resonances for the various contributions, is perhaps 
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the simplest to interpret [31, 32]. However, there are technical difficulties with broad 
and overlapping resonances in this approach. In principle the K-matrix formalism 
overcomes many of these difficulties [33], and analyses using it have begun to ap- 
pear [34]. However, there is relatively little experience with the K-matrix formalism 
in Dalitz analyses of charm decay. 

6.1.1 Dalitz Analysis of D+ — ^ tt+tt+tt" Decays from E791 




Figure 19: Dalitz plot projections of the E791 data and fits for tt+tt+tt" on the 

squares of invariant tttt masses: (a) for the fit without a low mass scalar contribution 
and (b) for the fit with a low mass scalar contribution. The histograms for the 
two ir'^TT' combinations are summed, and the shaded histograms are the background 
contributions. [10, 35] 



The E791 Dalitz analysis of tt+tt+tt^ decays utilizes 1,686 candidates pro- 

duced in 500 GeV/c 7r~-nucleon interactions [31]. In their isobar model they included 
contributions from the estabhshed resonances p°(770), /o(980), /2(1270), /o(1370), 
and p°(1450), and a non-resonant contribution. The projection of the fit with these 
resonances onto the 771"^+^- axis is illustrated in Figure 19(a). Clearly the ffi does not 
represent well the data in the low mass region, ^ 0.25 GeV/c^. Figure 19(b) 

shows that the fit is much closer to the data when E791 includes a low mass scalar 
cr" contribution. The fit fractions obtained by E791 are given in Table 12. It is 
interesting to note that the contribution has the largest fit fraction. 

6.1.2 Dalitz Analysis of £>+ tt+tt+tt" Decays from CLEO 

CLEO reported a preliminary Dalitz analysis of D'^ — >■ tt+tt^tt^ decay [36] using an 
isobar model similar to the one used by E791. The CLEO-c data sample included 
4,100 events with a signal to noise ratio of about two to one. Projections on the 
isobar ffi on the m'^+^- and 771^+^+ axes are illustrated in Figure 20. The preliminary 
CLEO results are also listed in Table 12. Except for the p(770)7r"'' contribution, the 
results of the two fits are in good agreement. 
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CLEO (%) 


E791 (%) 




41.8 ±2.9 


46.3 ±9.2 




20.0 ±2.5 


33.6 ±3.9 


/o(980)7r+ 


4.1 ±0.9 


6.2 ±1.4 


/2(1270)7r+ 


18.2 ±2.7 


19.4 ±2.5 


/o(1370)7r+ 


2.6 ± 1.9 


2.3 ± 1.7 


/o(1500)7r+ 


3.4 ± 1.3 




Non Res 


< 3.5 


7.8 ±6.6 


p(1450)7r+ 


< 2.4 


7.8 ±0.6 



Table 12: Fit fraction results of isobar Dalitz fits from E791 and CLEO for contribu- 
tions to -K^-K^-K" decays. The E791 results were published, while the CLEO 
results are preliminary. 

Z-projection 




All coherent 
X: Rho(770) 
X:f2(1270) 

X: Low mass pipi 



1 2 3 

m2(7i*7i-) (GeV/c^)^ 




m2(7i*7i*) (GeV/c^)^ 



Figure 20: Preliminary Dalitz plot projections on the squares of invariant vrvr masses of 
the CLEO-c data and fits for D'^ tt+tt+tt^: (left) projections on the vr+vr" invariant 
mass squared axis and (right) projections on the 7r+7r+ invariant mass squared axis. 
The empty bin at 'rn'^+^- = 0.25 GeV/c^ is due to a mass cut. 



6.1.3 Dalitz Analysis of D+ tt+tt+tt- Decays from FOCUS 

FOCUS reported [34] a Dalitz analysis performed on 1,527 tt+tt+tt" decay 

candidates. In the fit to the data, FOCUS included a nonresonant contribution, 
traditional isobar Breit-Wigner terms for p'^(770) and /2(1270) contributions, and an 
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Figure 21: Dalitz plot projections of the FOCUS data and fits for D"*" tt+tt+tt" on 
the squares of invariant tttt masses: (a) projections on the low-mass tt^tt" combination 
and (b) projections on the high-mass tt+vt" combination. [37] 



S-wave contribution. The S-wave contribution is not represented as a sum of scalar 
resonance contributions; it is represented by a K-matrix parameterization with five 
pole contributions. The projections of the high and low mass combinations of the 
FOCUS data and fits on the m^+^_ axis are illustrated in Figure 21. Perhaps the 
most illuminating comparisons of the FOCUS results with E791 and CLEO are the fit 
fractions of the p°(770) and /2(1270) contributions, which are 30.8±3.9 and 11.7±1.9, 
respectively. Although there is qualitative agreement among the three experiments, 
a comprehensive understanding of the resonant substructure of D"^ 7t~^7i~^7t~ will 
require substantially larger data samples and a better understanding of the strengths 
and weaknesses of the isobar and K-matrix formalisms. 

6.2 Dalitz Analyses that Contribute to Measuring 7 or 03 

Utilization of Dalitz analyses of decays [38, 39, 40] to enable measurements of 
the CKM angle 7 or 03 is one of the more surprising and interesting applications of 
these analyses. These techniques for measuring^ 7 utilize D'^(*^K^ decays, 

where Z)° can be or D and the superscript (*) denotes either D or D*. The 
core idea [41] is that the decay can be either Cabibbo favored [e.g., D^K^) or 

doubly Cabibbo suppressed {e.g., B~ D^K~). If the D'^ and decay to the same 
final state {e.g., — > i^^vr+vr") the amplitudes for these two decays can interfere. 
The phase of the interference term includes 7 and a strong phase; the former changes 

^To simplify notation, I will use 7 for either 7 or 03, or both. 
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sign between B~ and B'^ decays, while the latter does not. Hence, this interference 
can generate a direct CP violation, a difference between the rates of B'^ ID^^*^K^ 
and B~ D^^*^K~ decays. However, amplitude ratios and the strong phase in the 
D decay are required to determine 7 from a measured B^ decay rate asymmetry, and 
these parameters can be determined in a Dalitz analysis of or D decays to the 
appropriate final state. 

6.2.1 Dalitz Analyses oi D — > K^7z~^n~ Decays to Measure 7 

Recently, Babar [38] and Belle [39] analyzed the Dalitz plot for K^-k'^tt" decays 

in i?^ — s> D^^*^K^ decays. For the Dalitz analyses. Belle utilized 262,000 signal events 
from a 357 fb~^ data sample, and Babar utilized 81,500 signal events from 91.5 fb~^ 
of data. The purity of each final data samples was approximately 97%. Both 
groups fit their Dalitz plots with isobar models; BaBar fit 16 two-body states and 
Belle fit 18 two-body states, and the BaBar states are a subset of the Belle states. 
The projections of the data and fits on the Dalitz plot axes from the Belle and Babar 
are illustrated in Figures 23 and 22, respectively. Due to the enormous data samples, 
the errors on the data points are very small. The ability of the fits to match the 
data so well is remarkable. The fit fractions for the four most prominent resonant 
contributions are given in Table 13. Once again these analyses demonstrate that high 
quality charm decay measurements can be derived from the enormous BaBar and 
Belle data samples. 

From these analyses. Belle finds ^3 = 53° tigo ± 3° ± 9° and Babar finds 7 = 
70° ± 31° lio° -11°! where the third errors are estimates of the uncertainties in the 
Dalitz decay models. Within the large statistical errors there is good agreement 
between the two results. Since the statistical errors are larger than the systematic 
errors, there is room for improvement from the data that the two collaborations have 
not used in these analyses. 



State 


BaBar (%) 


Belle (%) 


JC(892)+7r- 


58.6 


61.2 




22.4 


21.6 




9.3 


9.8 


Non Res 


7.3 


9.7 



Table 13: The fit fractions for the four most prominent resonant contributions to 
— s> KgTT~^7r~ obtained by BaBar and Belle. 
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Figure 22: Dalitz plot projections on the squares of invariant masses of the BaBar 
data and fits for D — s> fiT^Tr+vr": (left) projections on the Kg7r~^ invariant mass 
squared axis, (middle) projections on the Kg7i~ invariant mass squared axis, and 
(right) projections on the vr+vr" invariant mass squared axis. [10, 42] 




ml {GeV^/C; 



(GeV%^) 



Figure 23: Dalitz plot projections on the squares of invariant masses of the Belle data 
and fits for D — >■ Kg7r~^7r~: (a) projections on the Kgir^ invariant mass squared axis, 
(b) projections on the Kg7r~ invariant mass squared axis, and (c) projections on the 
TT+TT" invariant mass squared axis. [10, 43] 



6.2.2 Dalitz Analysis of K+R-n^ Decays to Measure 7 

CLEO has studied the Dalitz plot of K+R-tt^ decays [40], which can also 

be used to measure 7 in —* D^K^ decays. The relative complex amplitude for 
K*^K~ and — ^ K*'^K~ decays is required to determine 7. This relative 
amplitude is the same as that for the two decays —>■ K*~K^ and ^ — *• K*^K~ 
(and their charge conjugates) assuming CP conservation in these decays. CLEO 
found 735 ^ K-R+tt^ candidates in 9.0 fb"^ of data taken with the CLEO III 
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detector. The Dalitz plot for the events were fit with 13 resonance components 
and a flat non resonant component. Four of these components, K*^ , K*~ , and 
non-resonant had the largest fit fractions. The projections of the fit on the three 
mass-squared axes are illustrated in Figure 24. The relative complex amplitude for 
the K*^K'^ decays is defined by, 

O'K'+K- 

where the a's are the real parts of the amplitudes and the J's are the phase shifts. 
The results for and 6d are = 0.52 ± 0.05 ± 0.04 and (f)D = 332° ± 8° ± 11°. 



3110406-002a 




0.6 1.0 1.4 1.8 0.6 1.0 1.4 1.8 1 2 3 



(mK+^o)' [GeV^/c^ (mK-7i°)' [GeV^/c^ (mK+K-)' [GeVVc^ 

Figure 24: Dalitz plot projections on the squares of invariant masses of the CLEO-c 
data and fits for —>■ K~^K~7r^: (a) projections on the K~^7r~ invariant mass squared 
axis, (b) projections on the K~ti^ invariant mass squared axis, and (c) projections on 
the K^K~ invariant mass squared axis. 

6.3 Analysis of K+R-tz+tz- Decays from FOCUS 

A study of resonant substructure in K^K~7r'^'K~ decay [16] from FOCUS is the 

only substantial effort to study resonances in four-body charm decays. After all cuts, 
FOCUS obtained 1279 ±48 K^K~t[^tt~ events above a modest background. A 

total of 10 resonant /decay contributions were considered in the analysis; three of these 
were A'i(1270)i^~ where the Ki decayed to states with other resonances, p°(770)-ft'"'", 
i<'Q(1430)7r+, and -R'*(890)7r+. The results of the fits projected onto invariant mass 
axes are illustrated in Figure 25. The fit fractions for the largest components are 
given in Table 14. 
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Figure 25: Two-body invariant mass distributions from the FOCUS Collaboration's 
study the four-body decay K^ti^tx^ . [44] 



Mode 


Fit Fraction (%) 




33 ± 6 ± 4 


i^i(1400)+ir- 


22 ± 3 ± 4 




29 ± 2 ± 1 


/s:*(1400)°i^+7r- 


11 ± 1 ± 1 


/o(980)7r+7r- 


15 ± 3 ± 2 



Table 14: Fit fractions for the largest components in the FOCUS Dalitz analysis of 
— s> K'^K^TT^TT" decays. 
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7 Summary and Conclusions 



The large data samples and high-quality detectors of the BaBar, Belle, and CLEO- 
c experiments have led to substantial advances in precision and discovery reach of 
charm physics. The CLEO Collaboration, operating the CLEO-c detector in the 
charm threshold region, is measuring absolute D hadronic branching fractions with 
unprecedented precision. 

CLEO reports preliminary results for and D'^ absolute branching fractions 
obtained from 281 pb~^ of e^e~ — *• '0(3770) — > D^D~ or D^D data. These results 
are limited by systematic errors that are as low as ^ 3%. CLEO expects to reduce 
these systematic errors via an ongoing effort. At this level of precision, final state 
radiation, whose effects are ^ 2%, must be considered. This is an interesting problem 
for the Particle Data Group, since most measurements of branching fractions do not 
take FSR into account. 

CLEO also reports preliminary results for hadronic Dg branching fractions from 
195 pb~^ of e'^e" annihilation data taken near Ecm = 4.17 GeV, which is just above 
the DfD*^ threshold. These results, with errors generally below 10%, are a substan- 
tial improvement over previous measurements, although they are limited by statistics. 
CLEO has an additional 130 pb~^ of data at this energy to be analyzed, and plans 
to take more data at this energy in the future. The resonant substructure of the 
Dg K~ K^TT^ decay mode is becoming an issue, since one of the contributions to 
this mode is 07r+, which has often been used as a reference branching faction 

for other decays. There appears to be a significant contribution (~ 5%) from 
scalar resonances in a reasonable region in the M{K^ K~) mass distribution around 
the peak. Since the error of the CLEO measurements of B{Ds — > K~ K^tt^) are 
comparable to the scalar contribution under the 0, we need to define a new reference 
branching fraction for Ds decays 

In addition to the measurements of Cabibbo-favored branching fractions men- 
tioned in the previous paragraphs, many new and accurate measurements of singly- 
and doubly- Cabibbo-suppressed branching fractions are appearing. BaBar, Belle, 
and CLEO-c are active in these analyses. By utilizing their enormous data samples, 
BaBar and Belle are starting to dominate measurements of the ratios of branching 
fractions to reference-mode branching fractions. 

There is substantially increased activity in the Dalitz analyses of hadronic final 
states of D decays. However, there are serious ambiguities in these analyses due to 
uncertainties in how to specify the amplitude for the resonant contributions. For ex- 
ample, E791, FOCUS, and CLEO have all reported Dalitz analyses of D'^ tt+tt+tt^. 
E791 and CLEO utilize a standard isobar parameterization of the decay amplitude, 
while FOCUS uses a K-matrix description which has some theoretical advantages. 
From these three experiments there is general agreement on the main resonant contri- 
butions to this decay, but there is important disagreement on the details, particularly 
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on the presence of a a scalar contribution to the decay. Developing a consensus on 
the proper treatment of resonant substructure with multiple overlapping states will 
become even more important as the larger data samples with higher purity become 
available. 

BaBar and Belle report Dalitz analyses of D K^-k'^ti" decays as byproducts of 
their studies of D^^*^K^ decays to measure the CKM angle 7 or 03. These 

collaborations are able to select enormous data samples with high purity from tagged 
B decays. CLEO reports a Dalitz analysis of K^K'-k^ decays which can also 

be used in a determination of 7 or 03. 

FOCUS reported the first analysis of the resonant substructure of a four-body 
hadronic charm decay, — > K^K'tt^-k^ . This result on a relatively modest data 
sample demonstrates that analyses of resonant substructure of charm decays to four 
or more bodies is a fertile field of research. 

In the near future, we can look forward to many more exciting and precise results 
from the charm sector from the BaBar, Belle, and CLEO collaborations. In the 
slightly farther future, we expect that the BESIII experiment operating in the charm 
threshold region at the BEPCII e^e~ storage ring in Beijing will open a new frontier 
in the precision and reach of charm decay experiments. 
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